Identification and Calculation of the Universal Maximum Drag Reduction Asymptote 

by Polymers in Wall Bounded Turbulence 
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Drag reduction by polymers in wall turbulence is bounded from above by a universal maximal drag 
reduction (MDR) velocity profile that is a log-law, estimated experimentally by Virk as V + (y + ) ~ 
11.7 log y + — 17. Here V + (y) and y + are the mean streamwise velocity and the distance from the wall 
in "wall" units. In this Letter we propose that this MDR profile is an edge solution of the Navier- 
Stokes equations (with an effective viscosity profile) beyond which no turbulent solutions exist. This 
insight rationalizes the universality of the MDR and provides a maximum principle which allows an 
ab-initio calculation of the parameters in this law without any viscoelastic experimental input. 



The mean streamwise velocity profile in Newtonian 
turbulent flows in channel geometries satisfies the clas- 
sic von-Karman "log-law of the wall" which is written in 
wall units as 



V + {y + ) = K~ 1 lny + + B, for y + > 30 



(1) 



Here x, y and z are the streamwise, wall-normal and 
spanwise directions respectively Q. The wall units are 
defined as follows: let p' be the fixed pressure gradients 
p' = —dp/dx, and L the mid-height of the channel. Then 
the Reynolds number Tie, the normalized distance from 
the wall y + and the normalized mean velocity V + (y + ) 
(which is in the x direction with a dependence on y only) 
are defined by 



Tie 



y 



(2) 



where vq is the kinematic viscosity. The law Q is univer- 
sal, independent of the nature of the Newtonian fluid; it 
is one of the shortcomings of the theory of wall-bounded 
turbulence that the von-Karmh constant ~ 0.436 and 
the intercept B « 6.13 are only known from experiments 
and simulations 0, Q . 

One of the most significant experimental findings Q 
concerning turbulent drag reduction by polymers is that 
in channel and pipe geometries the velocity profile (with 
polymers added to the Newtonian fluid) is bounded be- 
tween von-Karman's log-law and another log-law which 
describes the maximal possible velocity profile (Maxi- 
mum Drag Reduction, MDR) 0,0, HEf, 

V + (y+) = Kv - 1 ln(e« v2/ +) fory+>10. (3) 

This law, which had been discovered experimentally by 
Virk (and hence the notation k v ), is also claimed to be 
universal, independent of the Newtonian fluid and the 
nature of the polymer additive, including flexible and 
rigid polymers 8] . The numerical value of the coefficient 
k v is presently known only from experiments, k v _1 ~ 
11.7, giving a phenomenological MDR law in the form 

1 



V + (y + ) = 11.7 In y" 1 



17 



(4) 



For sufficiently high values of Tie and concentration of 
the polymer, the velocity profile in a channel is expected 
to follow the law J2J. For finite Tie , finite concentration 
and finite extension of the polymers one expects cross- 
overs back to a velocity profile parallel to the law 0J, 
but with a larger mean velocity (i.e. with a larger value 
of the intercept B) . The position of the cross-overs are 
not universal in the sense that they depend on the nature 
of the polymers and the flow conditions; the cross-overs 
are discussed in 0, • 

While we still cannot predict from first principles the 
parameters in von-Karman's log-law, the aim of this Let- 
ter is to identify the MDR log-law as an edge turbulent 
state in wall bounded flows, leading to a derivation of the 
parameters appearing in Eq. J3J without any viscoelas- 
tic input. The derivation follows from the theory of drag 
reduction by polymers that was developed recently, and 
therefore a short summary of the main aspects of the 
theory is in order. 

Wall bounded turbulence in Newtonian fluids is dis- 
cussed 0| by considering the fluid velocity U(r) as a 
sum of its average (over time) and a fluctuating part: 

U(r,t) = V(y)+u(r,t) , V(y) = (U(r,t)) . (5) 

The objects that enter the Newtonian theory are the 
mean shear S(y), the Reynolds stress W(y) and the ki- 
netic energy K{y) : 

S(y) = dV(y)/dy , W(y) = ~(u x u y ) , K(y) = (\u\ 2 )/2 . 

(6) 

In the presence of dilute polymers added to the New- 
tonian fluid one needs to complement these statistical 
objects with the dimensionless "conformation tensor" 
R(r,t) which stems from the ensemble average of the 
dyadic product of the end-to-end distance of the polymer 
chains, normalized by its equilibrium value Pq |lll Il2|. 
The way that the conformation tensor appears in the ad- 
ditional stress tensor which appears in the viscoelastic 
equations of motion is model dependent; it is different 
for flexible and rigid polymers, and it also depends on 
the actual model of the polymers. Nevertheless it was 



2 



shown theoretically 0, El that both for rigid and flexi- 
ble polymers one can write down eventually the balance 
equations for momentum and energy at distance y away 
from the wall as 



v{y)S{y) + W(y) = p'L , y^L, (7) 



K{y) = W(y)S(y) , (8) 
v{y) = v a + Cv p {R yv ) . (9) 



In Eq.Q the right hand side is the rate at which mo- 
mentum is generated by the pressure head, W(y) is the 
momentum flux in physical space towards the wall, and 
v{y)S{y) stands for the Newtonian viscous dissipation 
of momentum in addition to the polymer contribution 
to the dissipation of momentum. The effective viscosity 
u(y) is given by Eq. © where v v is the viscosity due to 
the polymers in the limit of zero shear and C is a con- 
stant of the order of unity. Similarly, in Eq. |jSJ the first 
term on the left hand side is the combined Newtonian 
and polymer contributions to the energy dissipation, the 
second models the inertial energy fluxes, and the right 
hand side is the (exact) energy production rate. The co- 
efficients a and b are dimensionless and of the order of 
unity. 

As in the Newtonian case, the balance equations need 
to be supplemented by a relation between K{y) and 
W(y). Rigorously the Cauchy- Schwartz ineqality leads 
to W(y) < K(y); experimentally one find that 



W(y) = c 2 v K(y) 



(10) 



with c v apparently y-independent outside the viscous 
boundary layer. To derive the functional form of the 
MDR [4( one asserts that the terms containing v{y) in 
the balance equations J7J and JHJ overwhelm the in- 
ertial terms, and then together with (|10() one derives 
S(y) ~ Const. /y which is the log-law for the MDR. Con- 
sisten with this law the effective viscosiy turns to be lin- 
ear in y. The increase in the viscosity of course increases 
the dissipation, but it was argued that the momentum 
flux W is reduced even further, leading to an increase 
in the mean momentum of the flow (for a given pressure 
head) , i.e. drag reduction. It is easy to argue 3 that the 
slope of the new log-law is larger that the slope in von- 
Karman's log-law, and hence drag reduction is obtained. 
Nevertheless, since neither c v nor the constants a and b 
in JSJ are known apriori, the actual slope of the MDR 
could not be determined. This shortcoming is remedied 
in the rest of this Letter. 

The crucial new insight that will explain the universal- 
ity of the MDR and furnish the basis for its calculation is 
that the MDR is a marginal flow state of wall-bounded 
turbulence. In other words this is the solution of Eqs. 
and © for which S(y) (or equivalcntly, the velocity 
profile) is the maximal possible for any turbulent solu- 
tion. Attempting to increase S(y)) further results in the 
collapse of the turbulent solutions in favor of a stable 
laminar solution W — 0. As such, the MDR is universal 



by definition, and the only question is whether a poly- 
meric (or other additive) can supply the particular effec- 
tive viscosity v(y) that drives Eqs. Q and ||SJ to attain 
the marginal solution that maximizes the velocity pro- 
file. We predict that the same marginal state will exist 
in numerical solutions of the Navier-Stokes equations fur- 
nished with a y-dependent viscosity u(y). There will be 
no turbulent solutions with velocity profiles higher than 
the MDR. 

To see this explicitly, we first rewrite the balance equa- 
tions in wall units. Define S +2 = a 2 K(y)/W(y), taken 
for simplicity as y-independent; we know from the New- 
tonian limit (in which v(y) = Vq) that S + ss 6 [l4j| . 
Once we write the equations with v(y) the ratio a 2 K/W 
can change drastically, and we denote it below by A 2 . 
With S+ = S v / p'L, K+ = K/p'L, W+ = W/p'L and 
u + (y + ) = v{y + )/v$. The balance equations are written 
as 



v + (y + )S + (y + ) + W+(y + ) = 1 
= S+ . 



A 2 

y + 



(11) 
(12) 



In Eq. A — » 6 + when v + {y + ) — ► 1 (the Newtonian 
limit) . The bunch of numerical constants in the second 
term on the LHS of l|12|l was replaced with k" 1 in agree- 
ment with newtonian log-law when v + (y + ) — > 1. In fact 
the second term on the left hand side of Eq. 112|) (which 
vanishes at the MDR) contains a factor (A/<5 + ) 3 . This 
factor is omitted for simplicity; accounting for this factor 
complicates slightly the algbra, leaving the final conclu- 
sions unchanged. Substituting now S + from Eq. I|ll|) 
into Eq. i(T2*)l leads to a quadratic equation for VW + . 
This equation has as a zero solution for W + (laminar so- 
lution) as long as v + (y + )A/y + = 1. Turbulent solutions 
are possible only when v + (y + )A/y + < 1. Thus at the 
edge of existence of turbulent solutions we find v + cx y + . 
This is not surprising, since it was observed already in 
previous work that the MDR solution is consistent with 
an effective viscosity which is asymptotically linear in y + 
0, 01 . It is therefore sufficient to seek the edge solution 
by maximizing the velocity profile with respect to linear 
viscosity profiles, and we rewrite Eqs. (| 1 1 11 and (|12|l with 
an effective viscosity that depends linearly on y + outside 
the boundary layer of thickness 5 + : 



[l + a{y+ -5 + )}S+ + W+ = 1 , 



[1 + a(y+ -£+)]- 



(13) 
(14) 



We now endow A with an explicit dependence on the 
slope of the effective viscosity, A (a) — > 6 + when a — > 0. 
For a / the ratio a 2 K/W is expected to depend on 
a (drag reduction involves a reduction in W), and this 
dependence is an important part of the following theory. 
We can present A (a) in terms of a dimensionless scaling 
function f{x), 



A(a) = S + f{aS+) 



(15) 
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FIG. 1: The solution for wVW+ (dashed line) and y + S + 
(solid line) in the asymptotic region y + ^> S + , as a function 
of a. The vertical solid line a = l/28 + = 1/12 which is the 
edge of turbulent solutions; Since VW+ changes sign here, 
to the right of this line there are only laminar states. The 
horizontal solid line indicates the highest attainable value of 
the slope of the MDR logarithmic law l/« v = 12. 



Obviously, /(0) = 1. In the Appendix we show that 
the balance equation (|13|) and (|14fl (with the prescribed 
form of the effective viscosity profile) have an non-trivial 
symmetry that leaves them invariant under rescaling of 
the wall units. This symmetry dictates the function A(a) 
in the form 



A(a) = 



6+ 



1 - a5+ 



(16) 



Armed with this knowledge we can now find the maximal 
possible velocity far away from the wall, y + 3> S + . There 
the balance equations simplify to 



ay + S + + W + = 1 



(17) 
(18) 



aA 2 (a) + VW+/k k = y+S+ . 

These equations have the y + -independent solution for 
v / W T and y+S+: 




+ 1 -a 2 A 2 (a) , 



y 



-S+ = aA 2 (a) + Vw+/n K 



(19) 



Obviously, (see Fig. the supremum of y + S + is 

obtained when W + vanishes, which happens precisely 
when a = 1/A(a). Using Eq. Ijlfijl we find the solu- 
tion a — a m = 1/2(5+. Then y + S + = A(a m ), giving 
k^ 1 = 25 + . Using the estimate 5 + w 6 we get the final 
prediction for the MDR. Using Eq. © with n^ 1 = 12, 
we get 



V + (y + ) w 121ny H 



17.8 



(20) 



This result is in close agreeement with the empirical law 
|J1| proposed by Virk. Note that the numbers appear- 
ing in Virk's law correspond to 5 + = 5.85, which is 



well within the error bar on the value of this Newto- 
nian parameter. Note that we can easily predict where 
the asymptotic law turns into the viscous layer upon the 
approach to the wall. We can consider an infinitesimal 
W + and solve Eqs. {H} and {TU for S + and the viscos- 
ity profile. The result, as before, is v + (y) = A(a m )y + - 
Since the effective viscosity cannot fall bellow the Newto- 
nian limit v + = 1 we see that the MDR cannot go below 
y + = A(a m ) = 25 + . We thus expect an extension of the 
viscous layer by a factor of 2. 

One should note that the result W + = should not 
be interpreted as W = 0. The differeence between the 
two objects is the factor of Tie 2 , W oc 7Ze 2 W + . Since 
the MDR is reached asymptotically as Tie — > oo, there is 
enough turbulence at this state to stretch the polymers 
to supply the needed effective viscosity. Nevertheless our 
discussion is in close correspondence with the experimen- 
tal remark by Virk that close to the MDR asymptote 
the flow appears laminar. 

In summary, one can probably improve further the 
model of the Newtonian wall-bounded flow, making it 
more elaborate and more precise. But the message of 
this Letter will stay unchanged; whatever is the model 
of choice, once endowed with effective viscosity v(y) in- 
stead of vq, there would exist a profile of v(y) that would 
result in a maximal possible velocity profile at the edge 
of existence of turbulent solutions. That profile is the 
prediction of the said model of choice for the MDR. In 
particular we offer a prediction for simulations: direct 
numerical simulations of the Navier-Stokes equations in 
a channel, endowed with a linear viscosity profile [5j, will 
not be able to support turbulent solutions when the slope 
of the viscosity profile exceeds the critical value that is 
in correspondence with the slope of the MDR. Notwith- 
standing, it is gratifying to discover that even a simple 
model of the balance of energy and momentum is suffi- 
cient, in light of the insight presented in this Letter, to 
predict ab-initio the functional form and the parameters 
that determine the Maximum Drag Reduction asymp- 
tote. 



APPENDIX A: THE SCALING FUNCTION 



Consider the following identity: 

v + {y + ) = l + a(y+ - S+) 



= [1 + a(y + - 6) + a(6 - <5+) 
a 



9(5) 



9(S) 



(Al) 



where 



g(S) = 1 + a(5 - S + ) , 6 > 5+ 



(A2) 
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Next introduce newly renormalized units using the effec- 
tive viscosity g(S), i.e. 



J = JL 



g(S) 



g{5) 



5* = — , S* = S+g{5) , W t = W+. 

(A3) 

In terms of these variables the balance equations are 
rewritten as 



[1 + a{y t - 8 t )}S i + W t = 1 
[1 + a{y t_ 5 t )] ^ 



y 



= s* 



(A4) 
(A5) 



These equations are isomorphic to ((13(1 and 1)14(1 with 5~* 
replaced by 8^ . The ansatz 1(15(1 is then replaced by 



Since the function A (a) cannot change due to the 
change of variables, the function A (a) given by Eq. (|A6(I 
should be identical to that given by Eq. (fT5|) : 



5+.f(a5+) = -^/K*) 



(A7) 



Using the explicit form of g(S) Eq. I)A2|) . and choosing 
(formally first) 5 = 5* = we find that /(£) = 1/(1 - £). 
It is easy to verify that this is indeed the solution of 
the above equation for any value of 8$ , and therefore the 
unique form of Eq. 1(161) . 



A(a) = -~rf(a8i) 
9{o) 



(A6) 



This form is dictated by the following considerations: (i) 
A(a) — * 8 + when a — > 0, (ii) all lengths scales in the 
rescaled units are divided by g (5), and thus the pre- factor 
in front of / becomes S + /g(S), and (iii) aS + in Eq. 1(12(1 
is now replaced in Eq. I(A5(I by aS$, leading to the new 
argument of /. 
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